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Представлен геотехнологический информационно-моделирующий комплекс, предназначенный для оптимизации добычи ура- 
на методом подземного выщелачивания. Комплекс включаете себя геотехнологическую информационную и моделирующую си- 
стемы. Геотехнологическая информационная система позволяет вводить, редактировать и отображать параметры, характери- 
зующие состояние продуктивного горизонта, а также характеристики технологических скважин. Моделирующая система описы- 
вает физико-химические процессы, происходящие в продуктивном горизонте при подземном выщелачивании урана. Комплекс 
может быть использован для оптимизации разработки месторождений урана методом подземного выщелачивания с целью по- 
вышения рентабельности разработки месторождений, увеличения доли извлекаемого урана и снижения экологической нагруз- 
ки на окружающую среду. 


Методом подземного выщелачивания (ПВ) раз- 
рабатывают экзогенные месторождения урана, ко- 
торые находятся в хорошо проницаемых подзем- 
ных водоносных горизонтах [1-3]. Извлечение 
урана из рудного тела происходит через систему 
технологических скважин. Через нагнетательные 
скважины в продуктивный горизонт нагнетается 
раствор веществ, способных растворять содержа- 
щие уран минералы. Образующийся в подземном 
водоносном горизонте продуктивный раствор из- 
влекается через откачные скважины. Образующие- 
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ся после переработки продуктивных растворов ма- 
точные растворы доукрепляются выщелачивающи- 
ми реагентами и снова подаются в нагнетательные 
скважины в качестве рабочих растворов. 

Основными задачами управления геотехнологи- 
ческим предприятием является повышение рента- 
бельности разработки месторождения, увеличение 
доли урана, извлекаемого из продуктивного гори- 
зонта, и снижение загрязнения подземных вод. Для 
решения этих задач нужно уметь оценивать остав- 
шиеся запасы урана, располагать информацией о 
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геохимическом состоянии продуктивного горизон- 
та и подземных вод, а также иметь возможность 
прогнозировать различные варианты развития 
предприятия и сравнивать различные способы раз- 
работки месторождения. То есть, необходима раз- 
вернутая во времени модель предприятия, позво- 
ляющая выбирать оптимальный способ разработки 
месторождения [4, 5]. В настоящей работе предла- 
гается геотехнологический информационно-моде- 
лирующий комплекс (ГТИМК), предназначенный 
для оптимизации добычи урана методом ПВ. 

В основу работы комплекса положено взаимо- 
действие геотехнологической информационной и 
моделирующей систем. Геотехнологическая ин- 
формационная система предназначена для отобра- 
жения пространственно-временных и атрибутив- 
ных данных на карте рассматриваемой территории. 
Визуальное представление информации базируется 
на цифровой модели продуктивного горизонта, ор- 
ганизованной в виде совокупности слоев и объек- 
тов, привязанных к единой системе координат. 
Под слоем понимается распределение какой-либо 
физической величины в рассматриваемой области. 
Объекты геотехнологической информационной 
системы представляют собой цифровые аналоги 
наблюдательных и технологических скважин и со- 
держат набор индивидуальных атрибутов-параме- 
тров, характеризующих реальный объект. Все слои 
и объекты геотехнологической информационной 
системы делятся на четыре модуля: геологический, 
гидрогеологический, минералого-геохимический и 
технологический. Геологический модуль содержит 
слои данных, характеризующих геологическое 
строение месторождения и физические свойства 
вмещающей породы (глубина залегания, мощность 
продуктивного горизонта, расположение рудного 
тела и т.д.). В минералого-геохимический модуль 
входят слои, включающие в себя данные о началь- 
ном распределении минералов в продуктивном го- 
ризонте и компонентов, растворенных в пластовых 
водах. Также в него включены кинетические коэф- 
фициенты и параметры равновесия, характеризую- 
щие физико-химические процессы, протекающие 
в системе. В гидрогеологический модуль включены 
слои, содержащие данные о распределении фильт- 
рационных параметров продуктивного горизонта 
(проницаемость среды, пористость, вязкость). 
Кроме того, в нем содержатся данные об условиях 
питания и разгрузки подземного водоносного го- 
ризонта. Объекты геотехнологической информа- 
ционной системы, являющиеся цифровыми анало- 
гами технологических объектов предприятия 
(скважин), сгруппированы в технологическом мо- 
дуле. Также в технологическом модуле содержится 
информация о режимах работы нагнетательных и 
откачных скважин и составах рабочих растворов. 

Моделирующая система основывается на ком- 
плексной математической модели многокомпо- 
нентной фильтрации и состоит из двух блоков, 
описывающих гидродинамические и физико-хи- 


мические процессы. Гидродинамический блок 
включает в себя расчет распределения давления и 
фильтрационных потоков в приближении жестко- 
го режима фильтрации. Гидродинамические расче- 
ты выполняются с учетом режимов работы техно- 
логических скважин, неоднородности фильтра- 
ционных параметров продуктивного горизонта и 
регионального потока подземных вод. В химиче- 
ский блок включены расчеты процессов растворе- 
ния урансодержащих минералов рабочим раство- 
ром, потерь кислоты в результате взаимодействия с 
кислотопоглощающими минералами, переотложе- 
ния минералов, комплексообразования, сорбции и 
десорбции. При расчетах учитывается неоднород- 
ность минералогического строения породы, кине- 
тика взаимодействия рабочего раствора с различ- 
ными минералами, неоднородность состава под- 
земных вод и переменный состав растворов, зака- 
чиваемых в продуктивный горизонт. 

Комплексная математическая модель описыва- 
ет фильтрацию многокомпонентного раствора в 
тонком ненаклонном пласте переменной мощно- 
сти (плановая фильтрация). Гидродисперсия учиты- 
вается с помощью приближения двойной пористо- 
сти [6]. Согласно данному подходу вся пористость 
представляется в виде суммы проточной /я, и за- 
стойной т 2 частей (т=т і +т 2 ). Проточными порами 
называются поры, по которым жидкая фаза способ- 
на двигаться под действием градиента давления. За- 
стойными порами называются поры, содержащие 
жидкость, не участвующую в адвективном движе- 
нии. Часть жидкости, находящаяся в проточных по- 
рах, называется мобильной, а в застойных - иммо- 
бильной. Рассматриваемая в модели система состо- 
ит из четырех частей: мобильная часть жидкой фа- 
зы (1), иммобильная часть жидкой фазы (2), вме- 
щающая порода, граничащая с мобильной частью 
жидкой фазы (3), вмещающая порода, граничащая 
с иммобильной частью жидкой фазы (4). Сжимае- 
мость жидкости и минерального скелета породы не 
учитывается, так как эти эффекты незначительны и 
не оказывают существенного влияния на абсолют- 
ную проницаемость пористой среды, а также не 
приводят к превращению проточных пор в застой- 
ные (или застойных пор в проточные). Скорость 
фильтрации 0 жидкости по проточным порам 
определяется законом Дарси [6, 7]: 

V = -— §гай (Р), 

№ 

где Р - давление; /и - вязкость жидкости; к - коэф- 
фициент проницаемости среды. Распределение да- 
вления рассчитывается с помощью закона нераз- 
рывности потока жидкости в приближении же- 
сткого режима фильтрации [8]: 

сііѵ(яс)=е, 

где Н- мощность пласта; 0 - плотность мощности 
источников и стоков жидкости. В качестве источ- 
ников и стоков жидкости выступают нагнетатель- 
ные и откачные скважины. 
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Процессы растворения и осаждения минералов 
в системе сернокислотный раствор - вмещающая 
порода, в общем случае, описываются уравнения- 
ми формальной кинетики для параллельных и по- 
следовательных химических реакций. Однако в на- 
стоящее время характер взаимодействия серной 
кислоты с большинством минералов, составляю- 
щих вмещающую породу подземных водоносных 
горизонтов, недостаточно хорошо изучен. Также 
неизвестны многие кинетические и равновесные 
параметры уравнений, описывающих эти взаимо- 
действия. В настоящей модели используется под- 
ход, основанный на выделении ограниченного чи- 
сла минералов, компонентов и процессов, которых 
достаточно для описания взаимодействия вмещаю- 
щей породы с серной кислотой. При этом реальные 
процессы, происходящие в подземных водоносных 
горизонтах, заменяются эффективными, описание 
которых позволяет учесть выход урана в раствор и 
убыль серной кислоты. В модели рассматриваются 
восемь групп минералов, которые различаются по 
составу (содержащие шестивалентный или четы- 
рехвалентный уран, кальций и др.) и скорости 
взаимодействия с раствором серной кислоты. При 
описании жидкой фазы рассматривается поведе- 
ние одиннадцати компонентов (1ГО 2 + , ІІ 4+ , Н + , 
80 4 2_ , Са 2+ , сульфатные комплексы уранила, урана и 
др.). Динамика концентраций всех компонентов в 
различных частях системы, за исключением иона 
водорода, описывается системой уравнений [9]: 
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где Сф - концентрация /-ого компонента в жидкой 
фазе (части системы Ф= 1,2); Щ- концентрация /- 
ого компонента, сорбированного на поверхности 
(Ф=3,4); СІ - содержание /-ого компонента в /-ом 
минерале в части Ф (Ф=3, 4); // 2 — плотность потока 
/-ого компонента из мобильной части жидкой фазы в 
иммобильную; /,/ (/ 2 ' 4 ) - плотность потока /-ого ком- 
понента из мобильной (иммобильной) части жидкой 
фазы на контактирующую с ней поверхность породы 
в результате сорбции (десорбции); (/$) - плот- 
ность потока /-ого компонента из мобильной (иммо- 
бильной) части жидкой фазы в у-ый минерал соот- 
ветствующей части твердой фазы в результате про- 
цессов растворения (/| ф+2 >0) или осаждения 
(/| Ф +2<0)- Концентрация иона водорода определяет- 
ся из условия электронейтральности раствора. Плот- 
ность потока // 2 между мобильной и иммобильной 
частями жидкой фазы принимается пропорциональ- 
ной разнице концентраций компонентов: 

2 12 — о: (С/ — С*2 ), 


где а 1 - постоянная скорости массообмена. Плот- 
ность потока /-ого компонента /| ф+2 из часть жид- 
кой фазы (Ф=1, 2) в граничащую с ней часть твер- 
дой фазы Ф+2 в следствие осаждения (растворе- 
ния) у-ого минерала определяется уравнением: 

4ф +2 = аіѲ(АГ ф )АГ ф +аі (ЯД 0(-. АР Ф )А Р ф (Ф=1,2) 

где АРф - движущая сила массопотока /-ого компо- 
нента; - постоянная скорости осаждения 

(растворения) у-ого минерала; (5 уд У Ф - удельная по- 
верхность у-ого минерала; Ѳ(х) - ступенчатая функ- 
ция (Ѳ(х)=0 при х<0 и Ѳ{х)=\ при х>0). В модели 
предполагается, что движущая сила А Р ф массопо- 
тока /-ого компонента в результате осаждения 
(АР>0) или растворения ( АР<0 ) минералов опре- 
деляется отклонением активности этого компо- 
нента в жидкой фазе от равновесного значения. 
Равновесные значения активностей компонентов 
являются функцией концентрации ионов водорода 
в соответствующей части жидкой фазы. Плотность 
массопотока /-ого компонента 4+2 в часть Ф жид- 
кой фазы из части Ф+2 твердой фазы в результате 
физической сорбции (десорбции) описывается в 
рамках приближения линейной кинетики: 

4 Ф +2 =-«И с 'р,ф ~С Ф ), (Ф=1, 2) 

где С/ф - равновесная концентрация і-ого компо- 
нента в растворе для процесса адсорбции (десорб- 
ции); а$ — константа скорости адсорбции (десорб- 
ции) /-ого компонента. 

ГТИМК создан на языке Вогіапсі Виіісіег С++ и 
представляет собой проблемно-ориентированное, 
многопоточное, многооконное, 32-битное програм- 
мное обеспечение, работающее под управлением 
операционной системы \Уіпс1о\у$ 98-ХР. Использова- 
ние средств интерфейса операционной системы \Ѵіп- 
бош при разработке ГТИМК предоставляет пользо- 
вателю возможность работы со стандартными объек- 
тами - окнами, списками, таблицами, которые легко 
связываются с базами данных и отображаются на 
экране монитора. Интерфейс ГТИМК включает в се- 
бя элементы управления - кнопки, переключатели, 
флажки, меню, всплывающие подсказки на русском 
языке. 

Запуск программы активизирует главное окно, 
содержащее заголовок, панель меню, панель инстру- 
ментов, рабочую область, статус-строку, отображаю- 
щую контекстные подсказки, расчетное время, коор- 
динаты положения курсора на карте, текущее си- 
стемное время. Главное окно приложения предназна- 
чено для работы с проектом, представляющим собой 
некоторую совокупность пространственно-атрибу- 
тивной информации о полигоне ПВ, геологическом, 
минералогическом строении и временной эволюции 
состояния продуктивного горизонта. Ввод и редакти- 
рование данных, управление расчетом и представле- 
ние результатов, сохранение данных осуществляются 
с помощью блоков управления проектом, управле- 
ния объектами геоинформационной системы и пара- 
метрами моделирования, визуализации данных, 
управления расчетом, обработки команд (рис. 1). 
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Блок управления объектами геоинформационной 
системы и параметрами моделирования предназна- 
чен для ввода и редактирования геологических, ги- 
дрогеохимических, минералогических и технологи- 
ческих данных о продуктивном горизонте, положе- 
нии и режимах работы скважин и составах закачива- 
емых растворов. Кроме этого возможно редактирова- 
ние параметров математической модели процесса 
ПВ. Блок управления расчетом позволяет запускать 
процесс моделирования, приостанавливать, остана- 
вливать или прерывать расчет, а также обеспечивает 
передачу данных в блок визуализации. Блок визуали- 
зации данных предназначен для отображения на кар- 
те местности скважин, реперов, сечений, изолиний, а 
также для управления режимом отображения данных 
в главном и дочерних окнах. Результаты моделирова- 
ния представляются блоком визуализации данных в 
виде картин пространственного распределения ис- 
ходных и текущих концентраций компонентов в 
жидкой и твердой фазах, содержаний минералов, па- 
раметров состояния продуктивного горизонта (пори- 
стость, мощность и др.) и профилей физических ве- 
личин вдоль задаваемых пользователем сечений. Воз- 
можен вывод временных зависимостей значений фи- 
зических величин в заданных точках (скважинах или 
реперах), а также технологических показателей по 
скважинам и по полигону в целом. 

Применение ГТИМК для проектирования и 
управления разработкой месторождения методом 


ПВ должно осуществляться в несколько этапов. На 
первом этапе создается цифровая модель продук- 
тивного горизонта. Предполагается, что к этому 
моменту проведены геологоразведочные работы, 
исследованы химические свойства вмещающей по- 
роды и подземных вод, накоплено достаточно ин- 
формации о состоянии месторождения. Имеющие- 
ся сведения, характеризующие минералогический 
состав продуктивного горизонта, его гидрологиче- 
ские и геологические свойства, а также другие па- 
раметры месторождения заносятся в базу данных. 
Далее генерируются распределения физических ве- 
личин во всей рассматриваемой области. Результа- 
ты экстраполяции анализируются специалистами, 
при необходимости производится пополнение ис- 
ходных данных и повторная генерация. На этом же 
этапе проводится определение параметров модели- 
рования физико-химических процессов. Параме- 
тры моделирования определяются по результатам 
лабораторных экспериментов и на основе данных 
опытных и опытно-промышленных работ. Резуль- 
татом первого этапа проектирования разработки 
месторождения должны являться цифровая модель 
продуктивного горизонта и набор физико-химиче- 
ских параметров моделирования. По результатам 
первого этапа подготавливается технико-экономи- 
ческое обоснование, содержащее результаты под- 
счета запасов урана, данные о фильтрационных и 
минералогических свойствах вмещающей породы 


81 



Известия Томского политехнического университета. 2005. Т. 308. № 3 


и выводы о целесообразности промышленного 
освоения месторождения. 

Если разработка месторождения признана рен- 
табельной, то нужно спланировать размещение 
технологических скважин. Создание цифровой мо- 
дели предприятия осуществляется на втором этапе. 
Цифровая модель предприятия включает в себя 
цифровую модель продуктивного горизонта и ха- 
рактеристики технологических объектов (нагнета- 
тельных и откачных скважин и т.д.). Наряду с пара- 
метрами технологических объектов в цифровую 
модель предприятия входят режимы их работы. 

На третьем этапе проводится прогнозный рас- 
чет эксплуатации месторождения. Основой для 
проведения расчета является цифровая модель 
предприятия. При расчете моделируются как ги- 
дродинамические, так и физико-химические про- 
цессы, протекающие в системе рабочий раствор - 
подземные воды - вмещающая порода. Результа- 
том расчета являются масса извлеченного урана и 
сопутствующих компонентов, расход реагентов 
(серная кислота и т.д.), оставшиеся запасы урана в 
продуктивном горизонте, а также другие характе- 
ристики разработки месторождения. На основании 
прогнозного расчета делаются выводы о целесооб- 
разности разработки месторождения данным спо- 
собом. Если способ не отвечает необходимым тре- 
бованиям, то проводится повторение третьего 
и/или второго этапов с целью поиска более прие- 
млемых схем разработки месторождения. На осно- 
ве совокупности проведенных прогнозных расче- 
тов подготавливается технико-экономическое 
обоснование и рабочий проект строительства гео- 
технологического предприятия, содержащие эко- 
номическую оценку разработки месторождения, 
прогнозные расчеты загрязнения подземных вод, 
предложения по проведению природоохранных 
мероприятий и другую информацию, необходимую 
для разработки месторождения. 

С помощью созданного геотехнологического ин- 
формационно-моделирующего комплекса было 


проведено моделирование разработки блока Дал- 
матского месторождения урана, Россия. На основе 
данных о минералогическом, геологическом и ги- 
дрогеологическом строении продуктивного гори- 
зонта в технологических скважинах, с помощью 
процедуры интерполяции, была создана цифровая 
модель блока месторождения. Физико-химические 
параметры определялись на основе сравнения ре- 
зультатов компьютерного моделирования и данных 
двухскважинноного опыта. Моделировалась разра- 
ботка блока, включающего в себя три откачных и 
восемь нагнетательных скважин. При этом исполь- 
зовались реальные характеристики скважин и режи- 
мы их работы. Распределение концентрации урана в 
жидкой фазе и водородного показателя рН в про- 
дуктивном горизонте, полученные спустя 1440 сут. 
после начала эксплуатации, приведены на рис. 2. На 
рис. 3 приведено сравнение интегральных техноло- 
гических показателей (массы откачанных урана и 
серной кислоты), полученных в ходе опытной раз- 
работки блока Далматского месторождения и ком- 
пьютерного моделирования. Хорошее соответствие 
расчетных и экспериментальных результатов под- 
тверждает адекватность модели и правильность рас- 
четов разработки месторождения методом ПВ. 

Разработан геотехнологический информацион- 
но-моделирующий комплекс, позволяющий моде- 
лировать разработку месторождения урана мето- 
дом подземного выщелачивания. Используемая в 
программном комплексе математическая модель 
адекватно описывает физико-химические процес- 
сы в продуктивном горизонте при сернокислотном 
выщелачивании урана. Программный комплекс 
может быть использован при проектировании и 
разработке месторождений урана методом ПВ для 
увеличения доли извлеченного урана, уменьшения 
расходов на единицу продукции и минимизации 
загрязнения подземных вод. На существующих ме- 
сторождениях, разрабатываемых с помощью ПВ, 
его можно использовать для обоснования ввода в 
эксплуатацию новых скважин, подбора режима эк- 



Рис. 2. Распределения концентрации урана в жидкой фазе (а) и водородного показателя рН (6) в продуктивном горизонте по- 
сле разработки блока Далматского месторождения в течение 1440 сут 
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Технические науки 




а) б) 

Рис. 3 . Сравнение показателей разработки блока Далматского месторождения и результатов компьютерного моделирования: 
а) масса извлеченного урана; б) масса откачанной серной кислоты, □ - экспериментальные данные, — - результаты 
моделирования 


сплуатации скважин, обеспечивающего макси- 
мальное извлечение урана при минимальных зат- 
ратах. На разведанных месторождениях ГТИМК 
может быть использован для определения опти- 
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